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ABSTRACT 

We find a simple, accurate model for the covariance matrix of the real-space cosmological matter 
power spectrum on slightly nonlinear scales (k ~ 0.1-0.8 /iMpc" 1 at z = 0), where off-diagonal matrix 
elements become substantial. The model includes a multiplicative, scale-independent modulation of 
the power spectrum. It has only one parameter, the variance (among realizations) of the variance 
of the nonlinear density field in cells, with little dependence on the cell size between 2-8 /i -1 Mpc. 
Furthermore, we find that this extra covariance can be modeled out by instead measuring the power 
spectrum of 8/a ce \\, i.e. the ratio of the overdensity to its dispersion in cells a few Mpc in size. Dividing 
6 by <7 ce ii essentially removes the non-Gaussian part of the covariance matrix, nearly diagonalizing it. 
Subject headings: cosmology: theory — cosmology: observations — large-scale structure of universe 
— methods: statistical 



1. INTRODUCTION 

The power spectrum of density fluctuations (and its 
Fourier dual, the correlation function) is the most com- 
mon statistic used to characterize the large-scale struc- 
ture of the Universe, and contains much information 
about cosmological parameters. In linear theory, the 
power spectrum's covariance matrix is diagonal, a prac- 
tical advantage over the correlation function. 

However, in the 'translinear' regime, the first decade or 
so in wavenumber smaller than fully linear scales, non- 
linearity grows in the density field. The dispersion cr ce n 
in the nonlinear overdensity field 5 in cells of translinear 
size 2ir/k ~ 4-60 h Mpc -1 increases from about 1 to 10 to 
20, as k decreases from 60 to 8 to 4 h Mpc -1 , as measured 
in the simulations analyzed below. (Note that the dis- 
persion in the linearly evolved density field in 8-/i _1 Mpc 
cells, cr 8 « 0.8.) 

The growing nonlinearity in 6 incites a substantial 
non-Gaussia n component in its power-spectrum covari 



ance matrix ([Meiksin 



1999 Takahashi eT~aT. 



White 



1999 
097 



Scoccimarro et al. 
1'he non-Gaussian 



covarian ce leads to a translinear plateau in Fisher infer 
mation ([Rimes fc Hamilton||200 5 2006 Neyrinck et al 



2006 Neyrinck &; Szapudi|2007p . When using the power 
spectrum of d to constrain a parameter such as the ini- 
tial power spectrum amplitude, modes in the translinear 
regime are highly correlated, giving little additional con- 
straining power when analyzed with larger-scale modes. 

A few methods have been proposed that mitigate this 
proble m, to varying degrees: remov ing large halos from a 



survey (Neyrinck & Szapudi 2007); applying a logarith- 
mic or rank-order-Gaussianizmg transform to the den- 



sity field before measuring the power spec trum (Neyrinck 
et al. 120091 |Seo et al.|2011l|Yu et al 120111); and nonlinear 



wavelet Wiener tittering (Zhang et al.||20ll[ ). 

The non-diagonality of the power-spectrum covariance 
matrix is not just a problem for parameter estimation. 
Without a theoretical model, off-diagonal elements must 
be estimated numerically, often cumbersomely. Even a 
covariance matrix estimated with 5000 simulations leaves 
residual noise (T09; see Fig. [4] below) . Also, the results 



often depend on the method used to measure the co- 
variance. For example, spatial jackknife r esampling sys- 
tema tically underpredicts the covariance ( jNorberg et al 
2009|. Another strategy, which we employ, applies a 



large set of sinusoidal weighti ng functions to a single den- 

HRS). 



sity field (Hamilton et al. 



2006 



The HRS weight- 
ings are designed to be optimally friendly to Fourier 
modes smaller than the sinusoids themselves, unlike win- 
dows with sharp configuration-space edges. Any window 
function introduces 'beat coupling' between small and 
large-scale modes; this contributes off-diagonal covari- 
ance beyond what is found in an ensemble of unweighted, 
periodic simulations. For the sinusoidal weightings, the 
beat coupling likely does not exceed that for more- 
disruptive, sharp window functions, often present in ob- 
servations. Another recent technique to generate low- 
noise covariance matrices from a relatively sm all number 
of simulatio ns makes use of 'mode-resampling' ( Schneider 
et al. 2011), taking advantage of the plethora of small- 
scale modes available. 

'Shrinkage estimation' (Pope & Szapudi 20081 is pos- 
sible if an approximate model is known for the covari- 
ance; this gives an optimal combination of a noiseless 
but idealized model, and a complete but noisy measured 
covariance matrix. One popular theoretical mode l for the 
covariance u ses a halo-model (HM) framework (Cooray 
& Hu |2001[ ). The full HM covariance includes up to 
tour-halo terms, and requires knowledge of the halo-halo 
power spectrum, bi- and trispectrum. Unfortunately, 
even the halo-halo power spectrum currently lacks a sat- 
isfying analytical model (Smith et al. 2007[ ). Leading- 
order perturbation-theory power, bi-, and trispectra can 
be used, but on translinear scales they are likely inac- 
curate for halo-halo correlations, making their nontrivial 
implementation possibly not worth the trouble. Thus, in 
practice often only the simpler 1-halo term is used. 

In this paper, we explore a simple, accurate model 
for transli near covaria n ce fro m a fluctuating multiplica- 
tive bias. Sato et al. (2009) considered a term with a 



rather similar form to this, called 'halo sampling vari- 
ance' (HSV), consisting of a variance in bias multiplying 
one-halo power-spectrum terms. They found that includ- 



2 






ing this HSV, which, unlike us, they ascribe to modes 
outside the survey volume, was necessary to model the 
covariance of the weak-lensing convergence power spec- 
trum accurately in a halo-model approach. 

2. RESULTS 

We measured real-space power spect ra and covariance 
matrices from the Coyote Univers e (Heitmann et al. 

^oT~ 



2010 
lations 



2009 Lawrence et al. 



Each of these has 



20101 ) suite of 37 simu- 
a different set of cosmo- 
each of them a plausible (given 



logical parameters 

current observational constraints) concordance ACDM 
model. The simulations occupy an orthogonal-array- 
Latin-hypercube in the five-dimensional parameter space 
n m h 2 6 [0.12,0.155], n b h 2 e [0.0215,0.0235], n s € 
[0.85,1.05], w e [-1.3,-0.7], a s e [0.61,0.9]. The 
remaining cosmological parameters, e.g. h, are set to 
match the tight cosmic microwave background (CMB) 
constraint on the ratio of the last-scattering-surface dis- 
tance to the sound-horizon scale. The 1024 3 -particle 
simulations have box size 1300 Mpc, fixed in Mpc (not 
h~ x Mpc) to roughly line up baryon- acoustic-oscillation 
(BAO) features in k among different cosmologies. Their 
resolution is sufficient for power-spectrum measurements 
accurate at sub-percent level at scales down to k = 
1/iMpc" 1 . We measured power spectra on grids with 
nearest-grid-point density assignment, not correcting for 
the pixel window function on the 128 3 -512 3 grids we use. 

The variable set of cosmologies prevented estimating 
covariances by directly comparing power spectra among 
simulations. So, we measured covariances within each 
simulation separately, using 248 (combinations of no 
weighting, and the first and second overtones in each 
Cartesian direction) sinusoidal weightings prescribed by 
HRS. We then averaged the covariance matrices from dif- 
ferent cosmologies together to reduce noise. Arguably, 
mixing cosmologies in this way is a good thing, decreas- 
ing the dependence on cosmological parameters. 

We did not use the covariance of the power spectrum 
P itself, but of InP, defining Cy = (AlnPjAlnP,) = 
(APiAPj) j (Pi) (Pj). Averaging this covariance ma- 
trix across cosmologies is more stable than the covari- 
ance matrix of P itself. We also define the quantity 
Tij = C l:j (N i N J ) 1 / 2 -26^, where 5% is a unit matrix, and 
Ni is the number of modes in bin i. Ty, zero for a Gaus- 
sian random field, is a measure of the non- Gaussian com- 
ponent of a covariance matrix, normalized to the Gaus- 
sian component. Conveniently for plotting purposes, it 
has roughly uniform noise properties. 
The top row of Fig. [I shows Ty, measured from 
1 1 the middle row shows, T|y is 

a = 



T tj = a{N t N fl 2 , for 



the simulations. As 
well-approximated by 

Var(cr^ ell )/cr^ ell . Here Var(cr^ cll ) is the variance of the 
measured (not linear-theory) cell-density variance cr? ell 
among sinusoidally weighted density fields, a, propor- 
tional to a variance, is inversely proportional to the vol- 
ume V = (1.3 Gpc) 3 /2 (the sinusoidal weightings effec- 
tively halve the volume). 

The approximate form of this covariance matrix is that 
of a Gaussian field, with additional power proportional 
to a scale-independent bias. Consider a field with power 
spectrum Pi — b (Pi) + P G . Here P G is the power spec- 
trum of a Gaussian random field such that (Pi) = (Pf), 



Cell size 
[Mpc/h] 



0.1 




Fig. 1. — Pixels are colored according to the quantities listed 
on the left edge of the plots; for quantities exceeding 1, contours 
appear at successive powers of 2. Top row. The normalized non- 
Gaussian part of the covariance matrix Ty. Middle row. An ap- 
proximation to Tij, using a = Var(cr? ,])/<7* |i, measured from each 
simulation's ensemble of sinusoidal weightings, and geometrically 
averaged across simulations. Bottom row. Tij for the power spec- 
trum normalized by the cell-density variance, P(k) / 'ciLji, equiva- 
lent to the power spectrum of (<5/<T ce n). (The displayed range of 
k includes modes not directly modulated by the weightings, and 
below the Nyquist k. The cell sizes are 1300/128, 1300/256 and 
1300/512 Mpc, roughly the listed values in h' 1 Mpc.) 

and b is a variable scale-independent bias, uncorrelated 
with fluctuations in P G , with (b) = 0. (Suppose b > — 1, 
e.g. is lognormally distributed, to avoid Pj < 0.) P's 
covariance matrix is 



C l3 = (b 2 )(P l )(P 3 ) + (APl 



; APf) 



(1) 



(APf APf) = 26(jl?/N u so = (b 2 ) (N^) 1 / 2 . Ne- 
glecting fluctuations from P G , cr^ ell = (b 2 ) a^ ell . Here 
the proportionality constant a = (b 2 ). 

With the insight that the covariance is from bias fluctu- 
ations, perhaps the covariance can be removed. Indeed, 
in the bottom row of Fig.[T] we show that the translinear 
covariance largely vanishes if the power spectrum from 
each density field is divided by cr^ e „, measured from that 
field. This is equivalent to measuring the power spectrum 
of (<5/c7 C eii). The reduction in covariance is dramatic for 
cell sizes > 4/i _1 Mpc. Using 2-hr 1 Mpc cells, there is 
significant residual covariance on small scales, but the 
covariance is still much smaller than in the top rows. In- 
terestingly, a is quite consistent over this factor of 4 in 
cell size; with 2, 4, and 8 ft." 1 Mpc cells, a x 10 4 =1.29, 
1.37 and 1.22. 

Fig.|2]shows a commonly plotted measure of the covari- 
ance, the correlation matrix CV,- / ' (CaCjj) 1 / 2 . Ps/cr^cn 
has smaller covariance even than the ln(l + S) power 
spectrum (NSS09), surpassed in diagonality only by the 
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Fig. 2. — Correlation matrices of various power spectra measured 
from the Coyote Universe simulations. Pg is the conventional power 
spectrum; Pg/o^n ls the power spectrum, dividing by the variance 
of cell densities; fwj+j) is the power spectrum of the log-density; 
PgCS) ls the power spectrum of the Gaussianized density. 



Fig. 4. 
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Fig. 3. — The non-Gaussian part of the covariance, , in the 
halo model, including just the 1-halo term (left panel), and with 
all terms (middle panel). as measured in the Coyote Universe 

simulations appears in the right panel. In all cases, Cg G is divided 
by a = Var(tr^ j,)/cr* ,j, measured from the simulations. 

Gaussianized power spectrum (NSS09). Hereafter, all of 
the Coyote Universe results use a 256 3 grid, with a cell 
size of ~ 4/i~ 1 Mpc. 

We note that the form of the ln(l+<$) covariance matrix 
is similar to the 8 covariance matrix, except that the con- 
stant a is given by fluctuations in the variance of ln(l+<5) 
in cells. This sheds light on the reduction in covariance 
from transforming the density to give a more-Gaussian 
1-point distribution. By definition, Gaussianization re- 
duces higher moments, in particular the kurtosis, and 
thus the variance in the variance. In our model, then, 
the full power-spectrum covariance is reduced. 

2.1. Comparison with the halo model 

Fig.gshows the quantity C^ G = C , 28% /(N.Nj) 1 / 2 , 
the covariance minus the Gaussian component; this 
should be a constant in our a model. We compare Cfj 
as measured from the Coyote Universe to predictions in 
the HM, both using all terms, and only the 1-halo term. 
For the HM covariance matrix we use our CosmoPy im- 



2008) 



plementatio n (|Neyrinck et al.||2006| |Neyrinck & Szapudi 
2007 



In all panels, 



and Tij measured from an ensemble of 5000 simulations by T09. 
The redshift-space Cf^ G is additionally multiplied by an empirically 

found factor <x (fc^fej) 1 / 2 to become roughly constant for large k. 
Note the linear k scale. 

by a as measured in the simulations. 

Indeed, at large fc, the measured C^ G /a w 1, conform- 
ing to the form above. For the HM panels, we set N% oc k 2 
exactly; i.e. fluctuations in Ni from the finite lattice are 
suppressed for clarity. The detailed features in the HM 
Cfj are different than that measured. On large scales, 

where there is much noise in the measured Cj? , there is 
barely a hint of t he perturbation-theory tris pectrum wig- 
gles from BAO flNeyrinck fc Szapudi|[2008| . This is not 
surprising; in that paper it took a tew hundred Gpc-scale 
simulations to get a clear signal. 

While our simple model involving a seems more accu- 
rate than the HM covariances on translinear scales, the 
HM does not fail utterly; the full-HM covariance is typ- 
ically within a factor of two of that measured. Given 
that we did not adjust HM parameters or polyspectrum 
assumptions to optimize the fit, its agreement is not 
bad. The agreement with the commonly used 1-halo co- 
variance, however, is poorer outside the large-fc corner. 
This vanilla 1-halo model assumes a halo mass function 
with uncorrclatcd Poisson fluctuations in each mass bin 



is normalized by dividing 



( Neyrinck et al. 2006 1 ; perhaps including mass- function 
covariance would improve agreement. 

2.2. Comparison with other results 

It is well worth checking our result with other covari- 
ance measurements, given our averaging procedure over 
cosmologies, and the 'beat-coupling' induced by our sinu- 
soidal weightings. Fig. [3] shows measurements by T09 of 
and c]] G from their suite of 5000 ACDM simulations 
of volume 1 (h^ 1 Gpc) 3 , each with 256 3 particles. 

For k > 0.04 ftMpc^ 1 , the real-space CfJ G is roughly 
constant, in accordance with our model. However, for 
smaller fc, C^ G roughly doubles from its smaller-scale 
value; unsurprisingly, other effects take over on the 
largest scales. As in the Coyote Universe results, at BAO 
scales there are interesting depressions in the covariance 
near the diagonal, where the covariance even dips slightly 
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The intrinsic Fisher information, (S/N) 2 



in various 



power spectra, measured from the Coyote Universe simulations. 
The dotted curve uses our model for the Pg covariance from fluc- 
tuations in a scale-independent multiplicative bias. 



negative. In Fig. [4] there is also a suggestion of periodic 
BAO bumps along the axes. 

In redshift space, we found empirically that C^ G in- 
stead falls off roughly as {kikj)^ 1 / 2 . We do not speculate 
here on why it has this form, but the reduced redshift- 
space covariance on small scales, noted by T09, is not 
surprising given the more-Gaussian redshift-space den- 
sity 1-point distributi on, smeared by redshift distortions 
(|Neyrinck et al.||2011[). 



2.3. Fisher information 

Fig.[5]shows signal-to-noise, (S/N) 2 , curves for the var- 
ious power spectra. This (S/N) 2 can be thought of as the 
effective number of independent modes over a range of 
wavenumber bins 1Z. 



(S/N) 2 (K) 



E (<#) 



(2) 



where C-n is the covariance matrix over 1Z. Here, the 
bins in 1Z vary from fc m i n , the largest modes not directly 
modulated by the sinusoidal weightings, to fc max . 

Compared to NSS09, the range of k plotted here is 
shifted smaller by a factor of 2 by the double box size, 
and there is reduced noise, but the trends are the same as 
found there. In fact, Ps/o- 2 cll is comparable to Pq(s), the 
power spectrum of the Gaussianized density. The (S/N) 2 
in Ps/&ce\i ram P s U P quickly at the smallest scales, be- 
cause the power spectra are effectively pinned together 
there when they are divided by c 2 eU , resulting in tiny 
(co) variance. We also show the (S/N) 2 curve for our 
Cfj G = a model, which approximates the Ps curve well. 

(S/N) 2 is a good single measure of intrinsic Fisher in- 
formation, but without derivative terms of the power 
spectrum with respect to cosmological parameters, the 
connection with parameter estimation is vague. For 
Ps, (S/N) 2 is roughly the information in In erf, un- 
marginalized over other parameters, since on linear scales 
d Ps (k) /(din a f) = 1, and it stays of order unity on non - 
linear scales, reaching ~ 2 (Neyrinck & Szapudi 2007 1 . 
However, at least for 8-/i -1 Mpc cells, Ps/cr'^ contains 



little information about erf. (If the measured power spec- 
trum were linear, obviously all sensitivity to cr| would be 
lost in dividing by it.) For Pg(<5); too, a free parameter is 
the variance in the Gaussianized density. In the present 
measurements, as in |Neyrinck et al. (2011 1, the G(S) vari- 
ance is fixed, pinning down the power spectrum at small 
scales in a similar manner as with P$/a£ on . The situation 



with Pi 



ln(l+<S) 



is in-between; on small scales it is allowed 



to fluctuate as it pleases (which is to a much lesser degree 
than Ps), but on large scales, < dPi n (i +( 5)/(<91no- 2 ) < 1 
because of large-scale bias pro duced by the log transfor m. 

An analysis of these issues ( Neyrinck|2011, in prep. ) is 
essential for these modified power spectra to be used for 
cosmological constraints. One might think that the in- 
sight that Pj's covariance is largely from fluctuations in 
a multiplicative bias might imply that this covariance is 
unimportant for parameters that depend on the power- 
spectrum shape, such as the tilt, or the BAO scale. How- 
ever, as shown in Fig. HI the form we have found for the 
translinear covariance does not really extend to fully lin- 
ear scales, leaving some room for shape fluctuations at 
the interface (occupied by BAO at z = 0) between the 
linear and translinear regimes. 

3. CONCLUSION 

We show that a fluctuating scale-independent multi- 
plicative bias provides quite an accurate model for the 
covariance matrix of the real-space dark-matter power 
spectrum on translinear scales. The non-Gaussian part 



of the covariance 



a 



Var(cr c 2 ell )/ t T4 ell , 



stant measure of the variance (among realizations) of the 
variance (within a realization) of the measured nonlinear 
density field, a is rather insensitive to the cell size used 
to measure it, at least within a range of 2-8 ft. -1 Mpc. a 
is not entirely trivial to measure, but it is much easier 
than the full covariance. For example, it could be esti- 
mated using sufficiently large simple sub-volumes, with 
little worry about edge effects. 

Furthermore, at least in real space, this translinear co- 
variance can largely be removed, or modeled out, by di- 
viding the power spectrum by cr 2 cll . This is equivalent to 
measuring the power spectrum of 5/a ce i\. However, this 
changes the sensitivity to cosmology, for example giving 
up information on the power-spectrum amplitude. 

The simplicity of this model suggests that it may be ap- 
plicable not only to dark matter, but perhaps to galaxies, 
or even to CMB anisotropics, where pushing to smaller 
scales dredges up some non-Gaussian covariance in the 
power spectrum o f Sunyaev-Zel'dovich effect anisotropics 
(Shaw et al.|2009). Shot noise would complicate the situ- 



ation for galaxies, but possibly this component could be 
accurately approximated with a Poisson term and sub- 
tracted. 
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